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ABSTRACT 

Using long-slit optical spectra obtained with the 2-m telescope at IUCAA Girawali 
Observatory, we show that the radio source J094221.98+062335.2 (z = 0.123) is as¬ 
sociated with a galaxy pair undergoing a major merger. Its companion is a normal 
star-forming galaxy infalling with a velocity of 185 km s _1 at a projected separation 
of 4.8 kpc. Using the Westerbork Synthesis Radio Telescope (WSRT) and Giant Me- 
trewave Radio Telescope (GMRT) we detect a strong H I 21-cm absorption at the 
systemic redshift of the radio galaxy with IV(H i)~ 9 x 10 21 cnU 2 for an assumed 
spin-temperature of 100 K. Such a strong H I 21-cm absorption is rare and has been 
seen only in a few compact radio sources associated with similar merging galaxy pairs. 
Milliarcsecond resolution Very Long Baseline Array (VLBA) observations resolve the 
radio source into a compact symmetric object with the hotspot separation of 89 pc. 
The 21-cm absorption is detected in the VLBA spectra towards both the radio lobes 
albeit with a strong optical depth gradient. We show that the strong 21-cm absorption 
is consistent with it being arising from a clumpy circumnuclear disk/torus. We also 
detect two weaker absorption lines redshifted with respect to the radio source in the 
WSRT/GMRT spectrum. They probably represent cold (i.e. T < 10 4 K) H I gas falling 
into the radio source. The presence of high concentration of H I gas in the circumnu¬ 
clear regions and signature of infalling cold gas allows us to conjecture that the young 
radio source may have been triggered by the gas infall due to the ongoing merger. 

Key words: quasars: active - quasars: absorption lines - quasars: individual: SDSS 
J094221.98+062335.2 - ISM: lines and bands 


1 INTRODUCTION 

It is now well recognised that active galactic nuclei (AGNs) 
play a vital ro le in the formation and evolution of their host 
galaxies (e.g . ISilk fc Reesl Il998l : iFerrarese fe Merrittl l2000l : 
iFabianl l2012h . The so-called AGN feedback through winds 
from the accretion disk or by the relativistic radio jets can 
deposit large amounts of energy and momentum into the 
interstellar medium (ISM) of the host galaxy and its envi¬ 
ronment, and is often invoked in theoretical and numerical 
simu lations to reproduce the observed properti es of galaxies 
(e.g. lCroton et al.1l2006l : ISchawinski et al.ll2006lb At galactic 
scales, radio jets can regulate or prevent star formation i.e. 
provide the most commonly observed negative feedback by 
expelling or heating the gas. The radio jets can also trigger 
star formation i.e. provide positive feedback by shocking the 
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gas to high densi ties thereby accelerating the gas cooling 
processes (see e. g. |yajmBrgugeLet_aL|[lh86; Chambers_etaL| 
Il987l : lDalvl[l990l : Reikuba et al. 20021 ; Gaibler et al. 2012f) . 


While all these observational results confirm that the 
AGN feedback strongly influences the evolution of its host 
galaxy and the environment, the conditions prevailing in 
the host galaxy may also quench or trigger the AGN ac¬ 
tivity by altering the rate of accretion onto the central su- 
permassive blackhole. One of the important questions in our 
understanding of the AGN is ‘What triggers the AGN activ¬ 
ity?’. The galaxy-galaxy mergers represent one of the natural 
mechanisms to funnel large quantities of gas to the central 
(<1 pc) regions of galaxies and trigger the AGN. Several au¬ 
thors have found evidence for connection between the merger 
process and the AGN activity for both r a dio-loud and 
radio-quiet AGNs (e.g. Combes et al. 2009; Ellison ^tjdJ 
201l]jRa : mqs_Almejda et_aJjj2£lj jViUnx-Ma^^ al.l 20121 : 


Tadhunter et al. 20121 ; Khabiboulline et al. 20141 '). However. 
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there are also studies that find no relationshi p be tween 
the AGN act ivity and mergers (e.g. ICisternas et al.l l201ll ; 
ISchmittll200ll 'l. In this context, it is of great interest to un¬ 
derstand in detail the conditions that trigger AGN in galaxy 
mergers. It has been suggested that cold gas associated with 
the circumnuclear disk or torus may play a key role in fu¬ 
eling the AGN activity. Therefore, a natural way to test 
the connection between galaxy-galaxy mergers and the trig¬ 
gering of AGN activity will be to observe the properties of 
circumnuclear gas in AGN hosts associated with mergers. 

For radio-loud AGNs, an interesting way of prob¬ 
ing the properties of circumnuclear atomic gas, i.e. its 
distribution and kinematics, is via H 1 21-cm absorp¬ 
tion observations of the compact steep-spectrum (CSS) 
and gigahertz peaked-spectrum (GPS) source^]. There 
have been several H 1 21-cm absorption searches in sam- 


sources fe.g. Carilli et al. 19981: 

van Gorkom et alJ 19891: 

Peck et all 200(4: iMor 

aanti et al. 

12001 

: Vermeulen et al. 

2003 

GuDta et al.l 2006 

Chandola et al. 

201 li: Allison et al. 

2012 

Chandola et al. 20131: iGereb et al.l 

2015h. High H 1 21- 


cm absorption detection rate (^30-50%), young age 
(<10 1 2 * * 5 yr) and the possibility of spatially resolved 21-cm ab¬ 
sorption to trace the parsec scale di stribution an d kinemat¬ 
ics of absorbing gas (see Table 1 of lArava et alJlioidl . for a 
nice summary) makes these sources excellent candidates for 
understanding the triggering of AGN and its early stages of 
evolution. 

As CSS and GPS sources are still embedded in their 
host galaxy ISM, they are also ideally suited for under¬ 
standing the impact of radio mode feedback in the form of 
radio jet-ISM int e raction fe.g. [Gelderman fe Whittlel 1994 


lAxon et al.l l2000l : iGupta et all 120051 : IShih et al.l 1201.4) , In 
the H 1 21-cm absorption studies, the AGN feedback due 
to CSS/GPS sources is revealed in the form of an ex¬ 
cess of blue-shifted absorption lines with outflow velocities 


sometimes as high as ~ 1000 km s 

Vermeulen et al.ll2003l; 

Morganti et al. 20051: GuDta et alJ 

20061: Chandola et al. 

201 ll: Mahonv et alj|2013; Morganti et al.ll2013i:lGereb et al. 

201511. The incidences of redshifted or infalling gas 


clouds fuelling the A GN are rare (e.g. lArava et al.l 20ld : 
iMaccagni et al.ll20l4) . 

In this paper, we report the detection and detailed anal¬ 
ysis of the H 1 21-cm absorption associated with a merging 
pair of galaxies at 2 ~ 0.123 having an angular separation 
of 2.2 arcsec. The absorption consists of a strong H 1 21-cm 
absorption component at the systemic redshift of the radio 
galaxy and two well-detached, infalling weaker components. 

This paper is structured as follows. In Section 2, we 
present details of optical and radio observations of the merg¬ 
ing galaxy pair serendipitously discovered through our on- 
going search for 21-cm a bsorption from quasar-galaxy pairs 
(IGupta et al . 1 I2010L l2013fPl In Section 3.1, we show that the 


1 The CSS and GPS objects typically have sizes of <15 and 1 kpc, 
respectively and believed to be you ng (<10 5 yr) counterparts of 
the larger radio galaxies and QSOs (lO'Poall 199iT i- 

2 Our sample of Quasar-Galaxy pairs at z < 0.3 is selected by 
cross-correlating galaxies detected in SDSS with radio loud AGNs 

in FIRST. There are 4 merging galaxy-galaxy pairs found. The 

21-cm and optical observations of the remaining 3 sources are 

underway. 


optical emission line ratios of the galaxy associated with 
the radio source, SDSS J094221.98+062335.2 (refer to as 
J0942+0623 from now on), are similar to those of AGNs 
whereas the optical spectrum of the other galaxy is consis¬ 
tent with a normal star forming galaxy. Besides a strong H 1 
21-cm absorption component at the systemic velocity of the 
radio-loud AGN, we also detect two redshifted absorption 
components representing a gas falling into the AGN. Sec¬ 
tion 3.2 focusses on deriving the properties of these H 1 ab¬ 
sorption clouds from the arcsecond and milli-arcsecond scale 
spectroscopy. In Section 3.3, we discuss the relationship be¬ 
tween the properties of H 1 21-cm absorption and merg¬ 
ers using a sample of compact radio sources with 21-cm 
absorption measurements. Finally in Section 4, we summa¬ 
rize our results. Throughout this paper we use the cosmo¬ 
logical model with Ha = 0.73, Qm = 0.27 and Ho = 71 
kms _1 Mpc _1 . 

2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Optical: IGO spectroscopy 

A single epoch optical spectrum of J094221.98+062335.2, 
observed on 27 February 2003, is available from the Sloan 
Digital Sky Survey (SDSS) archive. This spectrum covering 
a wavelength range of 3800-9200 A was obtained with a fibre 
of diameter 3” centered at the position of the radio source 
and contains flux from both the galaxies. Therefore, to de¬ 
termine the redshifts and the nature of both the objects, we 
performed the long-slit spectroscopy of J0942+0623 using 
the IUCAA Faint Object Spectrograph (IFOSC) mounted 
on the 2-m IUCAA Girawali Observatory (IGO) telescope. 
We used Grisms 7 and 8 of IFOSC0 in combination with a 
long slit of 1.5” width. These settings provide a wavelength 
coverage of 3800—6840 A and 5800—8350 A, and a spectral 
resolution of 300 km s _1 and 240 km s -1 , respectively. The 
Grism 8 data includes three 45 minutes exposures with the 
slit position angle, PA = 320° obtained on 25 February 2012 
and three 45 minutes exposures with PA = 0° obtained on 
26 April 2014. The Grism 7 data have one 45 minutes expo¬ 
sure with the PA = 0° obtained on 27 January 2011. The 
long slit orientations are shown in the left panel of Fig. [T] In 
the case of Grism 8, observations were carried out placing 
the object at different loc ations along the s lit to take care 
of the fringing effects fsee lVivek et al. 120091 . for details). 

For these data, the cleaning of raw frames and ID spec¬ 
tral extraction were carried out following standard proce¬ 
dures in irafQ. The two apertures that are shown overlaid 
on the 2-D spectra covering the H-a range in the top-right 
panel of Fig. [T] were chosen to extract ID spectra of the ob¬ 
jects, A (the host galaxy of the radio source) and B. The 
ID spectra covering the Ha and [N 11 ] lines are shown in the 
bottom-right panel. The spectra covering Ha, H/3, [N 11 ] and 
[O ill] lines overlaid with Gaussian fits are shown in Fig [2] 


3 Details of the IGO/IFOSC grisms can be found at 
http://www.iucaa ernet.in/~itp/etc/ETC/help.html#grism. 

4 IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1 . Left: SDSS colour composite image of the merging galaxy pair. The image, 20 ,/ X20 // in size, is centered at 0,12000 = 09h 
42 m 21.98 s, <5j2000 = 06 d 23 m 35.2 s. In the inset, we show the r-band image. The positions of two slits used for the IGO spectroscopic 
observations are marked. Top-right: The 2D spectrum covering the wavelength range of Ha emission for the exposures with slit PA = 
320°. The two apertures used to extract the ID spectra for the objects A and B are also marked. Bottom-right: One dimensional spectra 
corresponding to the two apertures showing the Ha and [N 11 ] emission lines. 


It is clear from these figures that the redshift of object A 
is less than that of B. We measure za = 0.1232 ± 0.0001 
and zb = 0.1239 ± 0.0001. The angular separation between 
the objects is 2.2 arcsec. This corresponds to a projected 
separation of 4.8 kpc at the redshift of the radio source. The 
redshift difference measured corresponds to a velocity sep¬ 
aration of ~185 km s _1 .The angular and velocity separa¬ 
tions found here are consistent with th e definition of merging 
galaxies used bv lEllison et alj (120111 1 to study the mergers 
in the SDSS database. In the following section, we study the 
line of sight distribution of the cold H 1 gas in this merging 
system using H 1 21-cm absorption line as a tracer. 

2.2 Radio: WSRT and GMRT spectroscopy 

We observed J0942+0623 with Westerbork Synthesis Radio 
Telescope (WSRT) on 3 March 2011 to search for the as¬ 
sociated 21-cm absorption. A total bandwidth of 10 MHz 
split into 2048 spectral channels was used to acquire the 
data in the dual polarization mode. This corresponds to 
a total bandwidth of ~2400km s -1 and a channel width 
of ~1.2 km s -1 prior to Hanning smoothing. The bright 
calibrators 3C 286 and 3C 147 were used for flux density 
and bandpass calibration. The total on-source time, ex¬ 
cluding calibration overheads, was ~9.8hr. The data were 
analysed using AIPi-0 following the standard procedures 
jGupta et al.ll20i(il f. We detect 21-cm absorption spread over 
~ 300 km s~ x in different components towards the radio 

5 See http://www.aips.nrao.edu/index.shtml for more informa¬ 
tion about AIPS 


source that is, compact in our WSRT maps having a res¬ 
olution of 180” x 14". The continuum normalised (stokes I) 
WSRT spectrum is shown in the top panel of Fig. [3] The 
measured continuum flux density is ~113mJy. The spectral 
resolution and the rms after applying Hanning smoothing 
are 2.4 km s -1 and 1.5 mjy beam -1 channel -1 , respectively. 
The velocity scale, defined with respect to the 2 em =0.1232 
corresponding to the systemic velocity of the radio source 
(i.e. object A), is also plotted on the top axis of the panel. 

Redshifted 21-cm absorption line data are often affected 
by radio frequency interference (RFI). Our WSRT observa¬ 
tions were no different in this regard. The data were care¬ 
fully examined to flag and remove the RFI. However, certain 
frequency ranges were affected by the interference that was 
present on almost all the WSRT baselines for most of the 
times. These frequency ranges are marked as shaded regions 
in the top panel of Fig.[3j and affect the detectability of weak 
absorption features in the range: 1260-1263 MHz. 

We observed the pair with Giant Metrewave Radio Tele¬ 
scope (GMRT) to confirm the weak H 1 absorption fea¬ 
tures detected in the WSRT spectrum. A bandwidth of 
4 MHz split into 512 frequency channels was used to cover 
the frequency ranges: 1259-1263 MHz and 1263-1267 MHz 
on 16 and 17 July 2011, respectively. In both the observing 
runs, 3C 147 was observed for flux and bandpass calibra¬ 
tion. In addition, nearby compact radio source 0925+003 
was observed every ~45 minutes for phase calibration. In 
each observing run the total on-source time was ~3hr. The 
data were analysed using AIPS following standard proce¬ 
dures. Relatively, the GMRT data were found to be much 
less affected by RFI. The flux density of the radio source 
J0942+0623 was found to be ~116mJy. Within the relative 
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Figure 2. IGO spectra of object A and B together with Gaussian fits to the H a , H 0 , [N II] and [O III] lines. It is clear that the 
emission lines are broader in the case of radio source (i.e. source A) which is at a lower redshift. 


uncertainties, this is consistent with the WSRT measure¬ 
ment. The radio source is unresolved in the GMRT maps 
with a resolution of ~3" x2". The stokes-I continuum nor¬ 
malised GMRT spectrum towards the radio source is shown 
in the bottom panel of Fig. [3] The spectral resolution is 
~1.9 km s -1 and the rms in the July 16 and 17 spectra are 
~1.7 and ~2.3mJy beam -1 channel -1 , respectively. Similar 
to the WSRT spectrum, the frequency ranges affected by 
RFI in most of the data are marked as shaded region in 

Fig- ED 

The velocity shift due to the heliocentric motion of 
Earth between the WSRT and GMRT observing runs is 
8 km s' 1 . This allows us to distinguish the real 21-cm ab¬ 
sorption components from the features due to RFI. A strong 
absorption component, labelled as ‘S’ in Fig. [3l and a weaker 
component, labelled as ‘Wi’, are detected in both the WSRT 
and GMRT spectrum at the observing frequency consis¬ 
tent with the heliocentric shift. Another component W 2 , al¬ 
though very close to the RFI ‘spike’ in the GMRT spectrum, 
is consistently reproduced in both the spectra. None of the 
other weaker features in the frequency range 1260-1263 MHz 


seen in the WSRT spectrum are reproduced in the GMRT 
spectrum and we consider these to be due to RFI. 

The properties of absorption components S, Wi and 
W 2 i.e. the absorption width, integrated 21-cm optical depth 
(f rdv) and H 1 column density (1V(H 1 )) are summarized 
in Table [T] Note that the inferred large JV(H l)~10 22 cm' 2 
is extremely rare in the samples of QSO sight lines trac¬ 
ing the ISM/halos of normal intervening galaxi es via 21-cm, 
Mg 1 1 or damped Lyman-q (DLA ) absorption (tGupt aetalJ 
1201(1120121 ; ISrianand et al.l 120121 INoterdaeme et al.l 120121 ') 


However, such N (H 1 ) a re more commonly ass ociated with 
Gamma-ray bursts (e.g. Ijakobsson et alJ |2006f), and galax¬ 


ies at an impact parameter of a few kpc (INoterdaeme et ahl 
l20l4) i.e. the regions of high density close to the cen¬ 
ter of galaxy. Amongst the known 21-cm absorbers asso¬ 
ciated with AGNs such column densities (i.e. ~10 22 cm' 2 ) 
have been seen in a few cases i.e. J0111+3906 (also 
known as OC314), SPS S J150805.7+342323.3, UGC6081 
UGC8387 and Mrk273 ( Carilli et al.l 199Sl : Chandola et al.l 
l201ll ; iDarling et all l201ll : Gereb et all 20151 1. It is clear 
from SPSS images that the host galaxies associated with 
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Figure 3. The WSRT (top) and GMRT (bottom) spectra of the radio source i.e. object A. The 21-cm absorption is clearly detected in 
a strong broad component (marked as S) at the systemic velocity of the radio source and two weak narrow absorption (marked as W 2 
and Wi) with an infalling velocity in the range ~80-200km s~ 1 . The velocity scale shown in the top panel is defined with respect to the 
systemic redshift 2 = 0.1232 of the radio source A. The vertical dashed lines correspond to redshifts of objects A and B. 


Table 1. Properties of the 21-cm absorption components as mea¬ 
sured from the WSRT spectrum 


Component 

Width’*’ 
(km s — *-) 

J r(v)dv 
(km s -1 ) 

A(H i)(100/T s )(/ c /1.0)t 
(10 20 cm" 2 ) 

s 

98 

49.9T0.4 

91.Oil.0 

Wi 

21 

0.54T0.09 

0.99i0.15 

w 2 

14 

0.63T0.07 

1.15i0.15 


t Velocity width corresponding to 90% of the total j rdv. 
t T s = spin temperature; f c = covering factor. 


J1508+3423 at 2 - 0.046, UGC6081 at 2 = 0.036, 
UGC8387 at z = 0.023 and Mrk273 at 2 = 0.037 are 
undergoing a major merger. In the case of higher redshift 
(z = 0.6685) source J0 111+3906, the host gala xy shows 
asymmetric morphology (jStanghellini et aDllOOSl ). possibly 
a signature of an ongoing or a recent merger. As galaxy 
mergers are expected to trigger the inflow of gas to the cen¬ 
tral regions of the galaxies, the large columns of H 1 gas 
associated with J0942+0623 and these five cases from the 
literature are caused by the ongoing mergers. 

2.3 Radio: VLBA imaging and spectroscopy 

The source J0942+0623 was observed with Very Long Base¬ 
line Array (VLBA) on 2012 June 15, 17 and 22. The goal of 
these observations was to perform the milliarcsecond resolu¬ 
tion 1264 MHz continuum imaging and the redshifted H 1 21- 


cm absorption spectroscopy. The total observing time was 
10 hr. The on-source time, excluding the telescope set-up 
time and calibration overheads, was ~6 hr. One 8 MHz- 
baseband channel pair covering the frequency range 1260.8- 
1268.8 MHz was used with right- and left-hand circular po¬ 
larizations sampled at 2 bits. The data were correlated with 
4000 spectral channels providing a spectral resolution of 
2 kHz i.e. ~0.47km s -1 . The correlator integration time 
was 2 s. We analysed the data following standa r d procedures 
in Al ps (see for example iMomiian et al.ll2002l : iGupta et ahl 
120121) . The continuum image of the radio source obtained 
by averaging absorption-free frequency channels is shown in 
Fig-El The rms in the image is 90 jrJybeam -1 , and the spa¬ 
tial resolution i.e. the restoring beam is 0.01230” x0.00436” 
with PA = —5.1°. 

In the VLBA image, we recover ~ 88% of the total flux 
density detected in the GMRT and WSRT images. Most of 
the milliarcsecond scale flux is contained in the two compo¬ 
nents which are marked with “+’ in the Fig. [4] Hereafter, 
we refer to these as the northern and southern blobs, or 
VLBA-N and VLBA-S, respectively. The separation between 
the peaks of VLBA-N and VLBA-S is ~40 mas i.e. 89 pc at 
z g =0.1232. The overall radio morphology is similar to the 
Compact Symmetric Objects (CS Os) believed t o be young 
(<10 5 yr old) radio sources fe.g.. IConwavl[2003 ). The peak 
and total flux densities of VLBA-N and VLBA-S are pro¬ 
vided in Tabled Both the blobs are spatially resolved and 
each can be fitted with two Gaussian components. We find 
that in each case more than 70% of the flux in the blob 
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Figure 4. VLBA image of J0942+0623. The rms in the image is 
90/Dy beam -1 . The contour levels are 4.Ox 10 -4 X(—1, 1, 2, 4, 8, 
...) Jy beam -1 . The gray scale range is: 0.4—30 mjy beam -1 . The 
restoring beam (see text for details) is shown as an ellipse at the 
bottom left corner. The image centre is at RA = 09 h 42 m 21 s .9832 
and Dec. = +06° 23 / 35” .1150. The northern and southern peaks 
are marked with ’+’. They are separated by ~40mas i.e. 89 pc. 


Table 2 . Results of VLB A observations 


Source 

Peak 

Integrated 

J rdv t 

component 

flux density 

flux density 



(mJy beam -1 ) 

( m Jy) 

(kms -1 ) 

VLBA-N 

31.0 

68.9 

33±3 

VLBA-S 

11.9 

31.4 

118±20 


f Towards VLBA-N and VLBA-S peaks. 


is contained in a Gaussian component coinciding with the 
peak with FWHM < 10 mas. 

In the top panel of Fig. [5] we plot VLBA spectrum ex¬ 
tracted towards the radio continuum detected in Fig. [4] The 
WSRT and integrated VLBA spectra match well within the 
uncertainties. This suggests that all the 21-cm absorption 
(except for Wi and W 2 which are too weak to be detected 
in the VLBA spectrum) detected in the WSRT/GMRT 
data are recovered in the VLBA data, and no absorption 
is present towards the radio flux resolved out in the VLBA 
data. The 21-cm absorption spectra towards the VLBA-N 
and VLBA-S peaks are shown in the bottom two panels 
of Fig. 0 As the optical depth sensitivity falls off rapidly 
beyond the continuum peaks, instead of making spatially 
resolved optical depth maps, we will make detailed compar¬ 
ison of the spectra towards the VLBA-N and VLBA-S peaks 
with our GMRT and WSRT spectra. These results are pre¬ 
sented in the next Section. 


3 RESULTS AND DISCUSSION 

In this section we present a detailed analysis of the optical 
and radio spectra and discuss the physical state of the merg¬ 
ing galaxies and the H 1 gas producing the 21-cm absorption. 


3.1 The merging galaxy pair: SDSS J0942+0623 


We use Gaussian fits to the emission lines to measure the line 
parameters towards the objects A and B. For the object A, 
we measure the FWHM of the Ha line to be 350 km s -1 . We 
find log [[O iii]A5007/H^]^0.4 and log [[N n]A6583/H a]~ 
—0.2. These values are consistent wi th the source A being 
a nar row line AGN (see Fig 14. 2 of lO sterbrock &; Fe rlandl 
l200Gl) . As per the definition of iKewlev et al.l (I2006I ). the 
source A can be classified as a composite H 11 -AGN type 
galaxy. In comparison, the emission lines associated with 
the object B are typically narrow and barely resolved in 
the IGO spectrum. We find log [[O ill]A500T/H/3]~ 0.0 and 
log []N n]A6583/H a] ~ —0.4. These values are consis¬ 
tent with what one expects from a H 11 region like object. 
Thus there are no signatures of any strong AGN activity 
in object B. Using the inferred Ha luminosity we estimate 
the dust uncorrected starfo rmation rate (using equations in 
lArgence fe Lamareilldl2009l ) of ~0.4 Mq yr 1 . This suggests 
that there is no strong star-bursting activities in the central 
regions of object B. Thus, the objects A and B form an in¬ 
teracting pair of galaxies, with a separation of 4.8 kpc, in 
which the former is a radio-loud narrow emission line AGN 
whereas the latter is a normal star-forming galaxy. 

It is not unusual for an AGN to be associated with 
a host galaxy undergoing a merger. Several authors have 
found evidence for mergers to be associat ed with both 
radio -l oud and rad i o-quie t AGNs (e.g. Comb es et al.l 
2009|_ Ellison et all 1201 ll : iRamos Almeida et al.l 20121 : 
VyLmMartm_e^al. 120121 : iTadhunter et al] 120121 : 

Khabiboulline et al.l 120141) . In a few merging pairs , 
double AGNs have als o been found (e.g. lEvans et al .1 l2008l ; 
ITadhunter et al.]|2012l ). Such systems are though extremely 
rare. In general, there could be several reasons for the 
lack of simultaneous AGN activities in both the merging 
galaxies like (i) intermittent supply of gas to the central 
regions or (ii) suppression of gas infall due to negative 
feedback caused by the nuclear starburst activity. In the 
case of object B, the lack of AGN activity is possibily due 
to the lack of cold gas supply to the central region. 


3.2 Properties of the absorbing clouds 

In this section, we first discuss the properties of component 
‘S’ inferred from the arcsecond and milliarcsecond scale 21- 
cm absorption spectra. Following this we present the prop¬ 
erties of well-detached narrow components Wi and W 2 that 
are infalling with respect to the systemic redshift of the ra¬ 
dio source. These features are detected only in the more 
sensitive GMRT and WSRT spectra. 


3.2.1 The 21-cm absorption component: S 

The 21-cm absorption component ‘S’ extends over the fre¬ 
quency range: 1264.4-1265.2 MHz. The absorption profile ex¬ 
hibits substructure indicating multiple absorption compo- 
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nents. As is evident in Fig [3] the absorption profile is asym¬ 
metric with respect to the systemic redshift of the radio 
source with the blue wing extending upto —120 km s' 1 and 
the red wing upto +50 km s' 1 . The 90% of the total op¬ 
tical depth is contained within ~98km s' 1 and the opti¬ 
cal depth weighted centroid of the profile is almost (within 
~20km s' 1 ) coincident with the systemic redshift of the 
radio source. 

Recently, iGereb et al.l (120151) used 21-cm absorption 
line shapes and systemic velocities for a sample of 32 ab¬ 
sorbers to study the origin of the absorbing gas. They sug¬ 
gested that while the narrowest lines with “busy function” 
FWHM<100 km s' 1 detected close to the systemic velocities 
are likely to be produced by H I disks, the absorption with 
intermediate and large (FWHM>200 km s' 1 ) widths origi¬ 
nate from H I disks with complex gas kinematics. Among the 
21-cm detections in their sample, three sources (UGC 05101, 
UGC8387 and Mrk273) are associated with mergers, and 
all three exhibit extremely broad absorption lines spreading 
over ~500-1000km s' 1 . Also, unlike J0942+0623 which is 
a CSO with an extent of 89 pc, the radio continuum asso¬ 
ciated with these mer gers have complex morphology with 
much larger extents ( Condon et al.l Il990l : ICrawford et al.l 
1 19961: ICarilli fc Tavlodl2000ll suggesting that the absorbing 
gas is originating from gas with complex kinematics due to 
ongoing merger and interaction with the radio source. 

The properties of component ‘S’ are akin to the class 
of H I absorbers with intermediate widths detected close to 
the systemic velocity. This suggests that the absorption orig¬ 
inates from a circumnuclear disk having simple kinematics. 
However, the narrower H I absorption line in our case can 
be attributed to the compactness of the J0942+0623 which 
may be tracing only a small portion of H I gas velocity dis¬ 
tribution. 

To understand this further, we compare the 
WSRT/GMRT spectra with the VLBA spectra towards the 
peaks of VLBA-N and VLBA-S in Fig. [o] The differences in 
the optical depth are clearly seen towards these two sight 
lines separated by 89 pc (see bottom two panels of Fig. 
The total integrated optical depth towards VLBA-S peak 
is higher by a factor of ~3.5 compared to that measured 
towards VLBA-N peak (see also Table [2]|. 

The top two panels in Fig. [o] show the GMRT and 
WSRT spectra. First we check whether the individual ab¬ 
sorbing gas detected towards the VLBA peaks cover all 
the radio emission from the corresponding VLBA compo¬ 
nent. For this we assume that the total flux seen in the 
WSRT/GMRT image is the sum of the fluxes of the two 
VLBA blobs. This allow us to predict the expected absorp¬ 
tion profile in the WSRT/GMRT spectrum i.e., 


= 0.314 exp{—rs{v)) + 0.686 exp(—T N (v)). 


(1) 


Here, I„ is the normalised flux measured with 
WSRT/GMRT and tn and rs are the optical depth 
measured towards VLBA-N and VLBA-S respectively. 
The values 0.314 and 0.686 are the fractional radio flux 
contained in the VLBA-S and VLBA-N respectively. 

The predicted absorption profile (solid smooth profile) 
using Eq. [lj is overlaid on the WSRT/GMRT spectra in the 
top two panels. It is clear that apart from the location where 
the peak absorption occurs towards VLBA-N, the predicted 
profile matches well with the observed one. This implies that 



1263 1264 1265 1266 


Observed frequency (MHz) 


Figure 5. VLBA spectrum towards the northern (VLBA-N) 
and southern (VLBA-S) peaks marked in Fig. [4] are shown in 
bottom two panels. Multiple Gaussian fit to the 21-cm absorption 
is also overlaid. In the bottom of each of these panel residual flux 
resulting from the fits are compared to the rms flux measured 
in the line free regions. Top two panels show the GMRT and 
WSRT spectra. In these, we overlay the total VLBA spectrum 
predicted using spectra towards VLBA-N and VLBA-S peaks (see 
text for details). The vertical dashed and dotted lines mark the 
emission redshifts and associated errors respectively of object A 
and B. The integrated VLBA spectrum extracted towards the 
radio continuum detected in Fig. 4 is shown as dotted line in the 
top panel. 


Table 3. Results of multiple gaussian fits to the 21-cm absorption 
detected towards the VLBA-N and VLBA-S peaks 


ID 

•^abs 

Tp 

FWHM 
(km s -1 ) 

rj^max 

kin 

(K) 

J rdv 
km s —1 


Fits for the profile towards VLBA-N 


Ni 

0.123245(1) 

0.98(0.10) 

17.04(1.27) 

6350 

17.8(0.3) 

n 2 

0.123159(2) 

0.29(0.03) 

30.22(9.02) 

21315 

19.6(0.5) 

n 3 

0.122975(5) 

0.24(0.04) 

19.70(3.59) 

8490 

5.0(0.2) 


Fits for the profile towards VLBA-S 


Si 

0.123154(1) 

1.72(0.16) 

9.36(1.33) 

1916 

17.2(0.5) 

S 2 

0.123139(2) 

1.56(0.05) 

50.36(1.70) 

55450 

83.6(0.6) 

s 3 

0.122940(4) 

0.56(0.03) 

26.38(2.13) 

15210 

15.4(0.2) 


the projected extent of the gas responsible for most of the 
absorption is at least as big as the individual continuum 
emitting regions in the VLBA images (i.e > 20pc). However, 
the above noted inconsistency over a narrow velocity range 
could be due to the presence of absorption by compact gas 
with dimensions smaller than 20 pc. 

Next we decompose the 21-cm H I absorption profiles 
towards VLBA-N and VLBA-S peaks using multiple Gaus- 
sians. We use these to examine if the properties of H I gas 
can be explained by a simple disk model. The fits are sum¬ 
marized in Table [3] and are plotted in the bottom two panels 
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Observed frequency (MHz) 


Figure 6 . Measured H I 21-cm absorption optical depth towards 
the two VLBA components are shown in the bottom two panels. 
Optical depth difference is shown in the second panel from the 
top. The top panel shows the fractional optical depth excess to¬ 
wards the southern component. Central frequency of individual 
Gaussian components are marked. The vertical dotted dashed 
lines are as in Fig [5] In the top panel the dotted horizontal lines 
give the mean fractional variation in r over the frequency range 
1264.6 MHz and 1265 MHz. 


of Fig. [ 6 ] Both the profiles are well fitted by three compo¬ 
nents. We label these components as Ni to N 3 and Si to S 3 
for VLBA-N and VLBA-S, respectively. 

It is apparent that a Gaussian component at z = 
0.123245 corresponding to Ni is not required to fit the profile 
towards VLBA-S. This probably corresponds to the narrow 
component over-predicted by the VLBA profile noted earlier 
and arises from a distinct cloud that at best covers only the 
VLBA-N. Similarly the component Si is seen only towards 
VLBA-S. For an assumed spin temperature of 100 K, the 
log N(H 1 ) in both the cases is ~21.5. If we assume the gas 
clouds to be spherical and covering only VLBA-N or VLBA- 
S then this will correspond to a gas density of nrr~50 cm' 3 . 
The density will be higher if Tsi§>100K as expected in the 
vicinity of a strong radio continuum source. For typical jet 
advancement speeds of ~ 0 . 1 c observed in the case of com¬ 
pact radio sources, the ambient medium with clouds of such 
densities can confine the radio source dTavlor et al.l l2000l: 


Bicknell et al.|[2003l : iGugliucci et al.ll200lj : Ijevakumar et all 
20051 ! . However, kinematically both Ni and Si are consis¬ 


tent with the systemic velocity of the host galaxy within 
~10km s' . This implies that these clouds are most likely 
associated with the circumnuclear disk or torus rather than 
the outflowing gas clouds interacting with the radio source. 

The other two Gaussian components fitted to the spec¬ 
trum towards VLBA-N and VLBA-S comprise ~75% of the 
total optical depth. These form a more coherent structure in 
opacity and kinematics that is consistent with a circumnu¬ 
clear disk/torus like structure. The components N 2 and N 3 


are within 10 km s ' 1 to the components S 2 and S 3 , respec¬ 
tively. This implies that there is no strong velocity gradient 
in the absorbing gas along the jet-axis i.e. the north-south 
direction. 

In top two panels of Fig. [ 6 j we plot the fractional optical 
depth, ( TS ~^ N ) an d the difference optical depth, (rs — tn) 
spectra. In the difference plot, we overlay the multiple Gaus¬ 
sian fit towards VLBA-S scaled by 0.85, the mean of ( TS ~^ N ) 
in the velocity range corresponding to S 2 and S 3 . Clearly 
the scaled rs fits match well with the difference spectra 
(i.e rs — tn)- This and the kinematical coincidence between 
N2 (N3) and S2 (S3) imply that these components cover 
both VLBA-N and VLBA-S albeit with an optical depth 
gradient from VLBA-N to VLBA-S i.e. over 89 pc. Follow- 
in g the circumnuclear disk model illustrated in the Fig. 5 
of lOrienti et al.l (l2006h . the opacity gradient across the ra¬ 
dio source can be explained by a combination of torus half¬ 
opening angle and, most importantly, the orientation of the 
line of sight such that the sight line towards VLBA-S tra- 
verses much longer path length through the d isk (see also 
iGupta fc Saikial 20061 : ICurran fc Whitingll2010l! . 

Thus, the VLBA spectroscopic observations are consis¬ 
tent with the absorption component ‘ S’ originating from a 
clumpy circumnuclear disk (as seen bv lMomiian et aDl2003l . 
in the case of NGC 7674). 


3.2.2 The H 1 21-cm. absorption components: W\ and W 2 

We detect 21-cm absorption component Wi and W 2 in the 
WSRT and GMRT spectrum. With respect to the radio 
source, these components have infall velocities of ~+85 and 
~+195km s ' 1 respectively. The widths corresponding to 
90% of the total optical depth and the integrated optical 
depths of these components are summarized in Table [T| The 
narrow velocity spread seen is consistent with the infalling 
gas being clumpy rather than a smooth infalling distribu¬ 
tion of gas where one expects to detect broad absorption 
lines. Typically, hotspot separ ation speeds of ~0 .1-0.5c have 
been observed for CSOs fe.g. iTavlor et aOl2000l l. Adopting 
a value of 0.1c for J0942+0623 we obtain the kinematic age 
of ~1500 yrs. This recent AGN activity may have been trig¬ 
gered by the infalling gas whose signatures are still seen in 
the form of Wi and W 2 . 

The component Wi at z a bs= 0.1235 is well repre¬ 
sented by a single Gaussian component of FWHM of 
6.4±1.0 km s' 1 . This places a constraint on the kinetic 
temperature to be less than 900 K and correspondingly 
log 1V(H l)^20. 910 . Similarly, the component W 2 at z a b a = 
0.12393 can also be fitted well with a single Gaussian compo¬ 
nent of FWHM of 11.5±1.5 km s' 1 . This places a constraint 
on the kinetic temperature to be less than 2900 K and cor¬ 
respondingly log N (H 1 )^ 21.46. From Tabic |T] it is clear 
that for any realistic spin-temperature the inferred N (H 1 ) 
will place the absorbing gas in the category of sub-DLAs or 
DLAs. The upper limit on the kinetic temperature derived 
for both these components are consistent with the infalling 
gas being cooler than what is expected (i.e. 10 4 K) for an 


6 Here we assume that the line width provides the maximum 
allowed temperature and spin temperature is equal to the gas 
kinetic temperature. 
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intergalactic H I gas in photoionization equilibrium with the 
UV background radiation. 

However, instead of being part of the infalling gas feed¬ 
ing the AGN, components Wi and W 2 could simply be from 
(1) an intervening tidally disrupted gas or (2) the gas associ¬ 
ated with objects A and B, respectively. This second possi¬ 
bility can be understood as following. As shown in Fig. [3] the 
blue wing of component ‘S’ extends upto +120 km s -1 . The 
infall velocity (+85 km s -1 ) of Wi is well within this range 
and it is possible that Wi also originates from the same 
circumnuclear disk as ‘S’. Similarly, the component W 2 ap¬ 
pears blueshifted with respect to object B by ~10km s . 
But this is within the errors (~30km s _1 ) on the redshift 
of object B implying that W 2 may just be associated with 
object B. 

3.3 Mergers and triggering of radio-loud AGN 

It has already been established through various 21-cm ab¬ 
sorption line surveys mentioned in Section 1 that detection 
rate of intrinsic 21-cm absorption amongst the compact ra¬ 
dio sources is significantly higher (~50%) than in the ex¬ 
tended sources (10-20%). This has been interpreted as the 
evidence of gas rich environment in these younger sources 
triggering the AGN activity. In this context, it is interest¬ 
ing of ask ‘what fraction of merging galaxies associated with 
compact radio sources are detected in H 1 21-cm absorption 
detections?’. 

This question can be addressed by con sidering the 
H 1 absorption study bv IChandola et ahl (1201 il l. The study 
is based on 18 CSS and GPS sources drawn from the 
CORA L Z sample ( Compact Radio Sources at Low Redshift; 
ISnellen et al.l 12004 1. As all the sources in this study are at 
low-z (z < 0.15) and 16/18 are also covered in SDSS, it is 
possible to determine the merger fraction for these sources 
by visual inspection. We find that three CSS/GPS sources, 
namely, ,11409+5216, J1508+3423 and J1718+5441 are as¬ 
sociated with galaxy-galaxy mergers and exhibit morpholog¬ 
ical distortions and tidal features typically showed by such 
systems. All three are detected in 21-cm absorption with 
H 1 column densities in the range: 10 21 ~ 22 cm~ 2 , which is 
at the higher end of column densities detected in their sam- 
ple. This is also tr ue for three mergers from the sample of 
iGereb et alj J2015l l discussed in Section 13.2.11 All this fur¬ 
ther reinforces the point already mentioned in Section 12.21 
that extremely high H 1 column densities (~10 22 cm~ 2 ), as 
seen in the case of ,10942+0623, appear to be associated with 
galaxies undergoing major mergers. It is important to note 
that, in general at 7V(H 1 )<10 22 cm -2 , one does see a mix of 
mergers and non-mergers, and the fraction of non-merging 
systems increases as the ff(H 1 ) cut-off is lowered. The 
major limitation at this point in quantifying this and draw¬ 
ing general conclusions is the unavailability of good quality 
optical imaging and spectroscopic data to reliably identify 
galaxy companions and mergers. 

The galaxy-galaxy merger timescales (~10 9 yrs) are 
substantially longer than the ages (<10 5 yrs) of radio sources 
considered here. Therefore, while mergers/interactions may 
be responsible for funneling large volumes of gas to the cen¬ 
ters of galaxies they may not be directly involved in trigger¬ 
ing the AGN activity. This is concurrent with the recent find¬ 
ings that while AGN fraction may be higher amongst merg¬ 


ing systems, mergers are not t he dominant cause of AGN 
activity fsee lEllison et allfeoill 'l. The AGN activity may ac¬ 
tually be triggered by local factors such as infall of gas due 
to disk instabilities. The 21-cm absorption components Wi 
and W 2 detected towards ,10942+0623 may be representing 
such gas clouds responsible for feeding the AGN. As noted 
in the previous section, the detection of infalling clumps in 
absorption is not sufficient to conclude that these are also 
responsible for fuelling the AGN. Additional constraints on 
gas geometry and kinematics, possibl y from emissio n lines 
from ionized and molecular gas (e.g. ICombes et al.l 120141 4 . 
are required to distinguish between various scenarios, and 
understand the role of mergers in triggering the AGN activ¬ 
ity. 


4 SUMMARY 

Using long-slit optical spectroscopy, with the 2-m telescope 
at IGO, we show that SDSS J094221.98+062335.2 is a merg¬ 
ing galaxy pair at z ~ 0.123 with a projected separation of 
4.8 kpc. One of these i.e. object A has radio emission typical 
of CSOs. The measured emission line widths and line flux 
ratios are consistent with this object being a narrow line 
AGN. The other galaxy, i.e. object B, is found to be a nor¬ 
mal star-forming galaxy and has an infall velocity of ~185 
kms -1 with respect to the radio source. 

We detect strong H 1 21-cm absorption associated with 
the radio source using WSRT and GMRT. The H 1 col¬ 
umn density, N (H l)=9.1x 10 21 (T S /100)(1.0// C ) cnU 2 is one 
of the highest amongst the associated 21-cm absorbers. 
We find that some of the strong 21-cm absorbers with 
N (H l)~10 22 cnU 2 reported in the literature are also asso¬ 
ciated with such major mergers. In the WSRT and GMRT 
spectra, we also detect two well-detached weaker 21-cm ab¬ 
sorption components that are redshifted with respect to ob¬ 
ject A by +85 and +195 km s _1 . These infalling absorptions 
could represent gas clouds fuelling the AGN. All this sug¬ 
gests a possibility that the large volumes of H 1 gas detected 
in these objects are caused by merger-driven inflows of gas 
that could probably be responsible for triggering and fuelling 
the central AGN. 

At milliarcsecond scale resolution, the radio source is 
resolved into a CSO like morphology with 88% of the arc- 
second scale flux in two dominant continuum components 
separated by 89 pc. The 21-cm absorption at velocities cor¬ 
responding to the strong 21-cm absorption component de¬ 
tected in the WSRT and GMRT spectra is detected towards 
both the continuum components albeit with an opacity gra¬ 
dient of as much as 3 over 89 pc. The VLBA spectra are 
not sensitive enough to detect the weaker redshifted absorp¬ 
tion components mentioned above. Comparing the VLBA 
with the WSRT/GMRT spectra we find that the properties 
of strong 21-cm absorption component are consistent with 
the absorption arising from a circumnuclear disk/torus as¬ 
sociated with object A. In this model, the opacity gradient 
across the radio source can be explained by a combination 
of torus half-opening angle and, most importantly, the ori¬ 
entation of the line of sight such that the sight line towards 
the southern VLBA component traverses much longer path 
length through the disk. 

The available data does not rule out the possibility that 
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instead of originating from gas inflow feeding the AGN, the 
absorptions at +85 and 195 km s _1 , could simply be from 
(1) intervening tidally disrupted gas or (2) gas associated 
with the ISMs of objects A and B, respectively. Additional 
constraints on gas geometry and kinematics are required to 
distinguish between various scenarios and understand the 
role of mergers in triggering the AGN activity. 
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